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ABSTRACT 

We have made a new calculation of the cosmic-ray secondary positron spectrum 
using a diffusive halo model for Galactic cosmic-ray propagation. The code computes 
self-consistently the spectra of primary and secondary nucleons, primary electrons, 
and secondary positrons and electrons. The models are first adjusted to agree with 
the observed cosmic-ray Boron/Carbon ratio, and the interstellar proton and Helium 
spectra are then computed; these spectra are used to obtain the source function for 
the secondary positrons/electrons which are finally propagated with the same model 
parameters. The primary electron spectrum is evaluated, again using the same model. 
Fragmentation and energy losses are computed using realistic distributions for the 
interstellar gas and radiation fields, and diffusive reacceleration is also incorporated. 
Our study includes a critical re-evaluation of the secondary decay calculation for 
positrons. 

The predicted positron fraction is in good agreement with the measurements up 
to 10 GeV, beyond which the observed flux is higher than that calculated. Since the 
positron fraction is now accurately measured in the 1-10 GeV range our primary 
electron spectrum should be a good estimate of the true interstellar spectrum in this 
range, of interest for gamma ray and solar modulation studies. We further show that 
a harder interstellar nucleon spectrum, similar to that suggested to explain EGRET 
diffuse Galactic gamma ray observations above 1 GeV, can reproduce the positron 
observations above 10 GeV without requiring a primary positron component. 



Subject headings: cosmic rays — diffusion — elementary particles — Galaxy: general 
— ISM: abundances — ISM: general 



1. Introduction 

Secondary positrons in galactic cosmic rays are an important diagnostic for models of 
cosmic-ray propagation and also for solar modulation studies. Recently several new experiments 
have provided improved data on both the positron to electron ratio and the positron spectrum 
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itself ( [Barwick et al. 1997 , Barbiellini et al. 1996| , Golden et al. 1996). These and previous 



data have mainly been compared with model predictions using leaky-box and diffusion models 
from Protheroe (1982). Since these predictions were made, new information and ideas have 
contributed to our understanding of cosmic-ray propagation: new measurements of nucleon 
secondary-to-primary ratios, the energy dependence of the propagation, results from gamma-ray 
astronomy and studies of the effects of diffusive reacceleration. For this reason new calculations of 
secondary positrons are desirable. 

Recently an extensive computer code for the calculation of galactic cosmic-ray propagation 
has been developed ( [Strong fc Moskalenko 1997a| ), which is a further development of the 
approach described by Strong & Youssefi (1995) and Strong (1996). Primary and secondary 
nucleons, primary and secondary electrons, and secondary positrons are included. The basic 
spatial propagation mechanisms are (momentum-dependent) diffusion and convection, while in 
momentum space energy loss and diffusive reacceleration are treated. Fragmentation and energy 
losses are computed using realistic distributions for the interstellar gas and radiation fields. 

The main motivation for developing this code was the prediction of diffuse Galactic gamma 
rays for comparison with data from the Gamma Ray Observatory CGRO instruments EGRET, 
COMPTEL and OSSE. More generally the idea is to develop a model which self-consistently 
reproduces observational data of many kinds related to cosmic-ray origin and propagation: directly 
via measurements of nuclei, electrons and positrons, indirectly via gamma rays and synchrotron 
radiation. These data provide many independent constraints on any model and our approach 
is able to take advantage of this since it must be consistent with all types of observation. We 
emphasize also the use of realistic astrophysical input (e.g. for the gas distribution) as well as 
theoretical developments (e.g. reacceleration). The code is sufficiently general that new physical 
effects can be introduced as required. In this paper we focus on positrons. The secondary positron 
problem is closely connected with that of gamma-rays from 7r°-decay since the same cosmic-ray 
nucleons are involved. For this reason the present application is a natural extension of the original 
code. 

The basic procedure is first to obtain a set of propagation parameters which reproduce the 
cosmic-ray B/C ratio and the spectrum of primary electrons. The proton and Helium spectra are 
then computed using these parameters, and used to obtain the source function for the secondary 
positrons and electrons. The secondary positron and electron spectra are then computed using 
the same propagation model. Note that although the ratio of positrons to electrons is normally 
used for the comparison with data, it is in fact the positron spectrum itself which is of interest in 
testing propagation models; however since the ratio is more easily measured experimentally we 
also use it in this work. 

Turning the argument around, we note that if we can compute the positron spectrum 
reliably, then accurate measurements of the positron fraction enable us to determine the true 
interstellar electron spectrum unless there is significant charge-sign dependent solar modulation 
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(for a discussion of the charge-sign dependence see Clem et al. 1996 ). This approach also provides 
an independent test for models of the primary electron spectrum, complementing gamma-ray 
and synchrotron measurements. Especially around 1 GeV the positron fraction is now rather 
accurately known, so the interstellar electron flux at this energy can be determined. 

We note that similar calculations of the electron spectrum have been made by Porter & 
Protheroe (1997) using a (1-D) Monte Carlo approach. They also address the galactic diffuse 
gamma-ray spectrum. However, since the main subject of the present paper is positrons, a detailed 



comparison of our electron and gamma-ray spectra with theirs will be made elsewhere ( Strong fc 



Moskalenko 1997b). 



2. Description of the models 



The models will be described in full detail elsewhere ( Strong fc Moskalenko 1997b| ; see also 



the information in Section 7); here we summarize briefly their basic features. The models are three 
dimensional with cylindrical symmetry in the Galaxy, and the basic coordinates are (R, z,p) where 
R is Galactocentric radius, z is the distance from the Galactic plane and p is the total particle 
momentum. The distance from the Sun to the Galactic centre is taken as 8.5 kpc. In the models 
the propagation region is bounded by R = Rh, z = Zh beyond which free escape is assumed. We 
take Rh = 30 kpc. The case Zh = 3 kpc has been studied since this is consistent with studies of 



radioactive nuclei (the Be/Be ratio: Lukasiak et al. 1994]) and synchrotron radiation. For a given 



Zh the diffusion coefficient as a function of momentum is determined by B/C for the case of no 
reacceleration; with reacceleration on the other hand it is the reacceleration strength (related to 
the Alfven speed) which is determined by B/C. Reacceleration provides a natural mechanism to 
reproduce the B/C ratio without an ad-hoc form for the diffusion coefficient ( |Simon fc Heinbach 



1996 , Heinbach fc Simon 1995 , 5eo fc Ptuskin 1994 , Letaw, Silberberg, fc Tsao 1993|) . The spatial 



diffusion coefficient for no reacceleration is taken f3D$ below rigidity po, f3D$(p/ po) s above rigidity 
Pq. The spatial diffusion coefficient with reacceleration assumes a Kolmogorov spectrum of weak 
MHD turbulence so D = (3Dq(p/ po) s with 5 = 1/3 for all rigidities. For the case of reacceleration 
the momentum-space diffusion coefficient D pp is related to the spatial coefficient using the formula 
given by Seo and Ptuskin (1994) (their equation [9]), and Berezinskii et al. (1990). The main free 
parameter in this relation is the Alfven speed va- The injection spectrum of nucleons is assumed 
to be a power law in momentum for the different species, dq(p)/dp tx p -7 for the injected density, 
corresponding to an injected flux dF{p)/dp oc f3p~" f or dF(Ek)/dEk tx p~ y since dp/dEk = A/f3, 
where E^ is the kinetic energy per nucleon, (3 = v/c. The value of 7 can vary with species. The 
interstellar hydrogen distribution uses information from HI surveys (|Gordon fc Burton 1976j , |Cox 



Kriigel, fc Mezger 1986j ) for the atomic component and CO surveys ( |Bronfmann et al. 1988j ) for 



the molecular hydrogen; we use also information on the ionized component. The Helium fraction 
of the gas is taken as 0.11 by number. The interstellar radiation field for inverse Compton losses 
is based on stellar population models and IRAS ( Bloemen, Deul, fc Thaddeus 1990| ) and COBE 
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(Boulanger & Desert 1992) analyses; in addition we include the cosmic microwave background. 
The magnetic field is assumed to have the form Bj_ = 6 e -\ z \/ 5k P c - R / 20k P c ^Q, This magnetic 
field model is based on recent observational work (Vallee 1994, Heiles 1996j ), and our local field 



of 4 jiG is consistent with these papers. The z-variation was chosen to follow studies of radio 



synchrotron surveys, which suggest a vertical extent of emissivity of several kpc (Phillipps et al 



1981 ) and hence a comparable extent of both field and electrons. Energy losses for electrons by 
ionization, Coulomb interactions, bremsstrahlung, inverse Compton and synchrotron are included. 
The inverse Compton losses were computed using the Thomson approximation, which is sufficient 
for our purposes. The Klein-Nishina corrections affect only electrons of energies above 100 GeV 
interacting with optical photon fields, while we are interested here in electron and positron energies 
below 100 GeV, and in any case far-infrared photons dominate the interstellar radiation energy 
density. Energy losses for nucleons by ionization and Coulomb interactions are included following 
Mannheim & Schlickeiser (1994). The distribution of cosmic-ray sources is chosen to reproduce 
the cosmic-ray distribution determined by analysis of EGRET gamma-ray data ( |Strong &: Mattox 



1996). 



First, the primary propagation is computed giving the primary distribution as a function of 
(R,z,p); then the secondary source function is obtained from the gas density and cross-sections, 
and finally the secondary propagation is computed. The entire calculation is performed with 
momentum as the kinematic variable, since this greatly facilitates the inclusion of reacceleration. 
Nucleon spectra are afterwards converted to dF(Ek)/dEk for comparison with observations. 



3. The model parameters 

Table |l] lists the parameters which were used for the models without and with reacceleration. 
The cross-sections for secondary production from the progenitor C, N, O were taken from Webber, 
Lee, & Gupta (1992), Heinbach & Simon (1995) and references therein. The source relative 
abundances were taken from Engelmann et al. (1990). The injection spectrum for Carbon was 
taken as dq(p)/dp oc j> -2,35 , for the case of no reacceleration, and p~ 2 - 25 with reacceleration. These 
values are consistent with Engelmann et al. (1990) who give an injection index 2.23 ± 0.05. 

Figure [I] shows the predicted and observed B/C ratio for the adopted parameters. We use 
the Voyager data from Webber et al. (1996). The spectra were modulated to 500 MV appropriate 
to the Voyager data using the force- field approximation flGleeson fc Axford 1968| ). The agreement 
is sufficiently good for our purpose of computing the secondary positrons. We can note that 
below 1 GeV the better fit is given by the reacceleration model (which also has one less free 
parameter), but this has no effect on the positron calculation. Inclusion of a convective term 
would be necessary to improve the fit in this region for the no reacceleration case. 

The proton and Helium spectra are computed as a function of (R,z,p) by the propagation 
code. Fig. |2| shows the predicted and observed proton and Helium spectra. The injection spectrum 
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is adjusted to give a good fit to the locally measured spectrum, normalizing at 10 GeV/nucleon. 
For protons we use Seo et al. (1991: Fig. 10a) based on LEAP and IMP8 balloon measurements, 
and Mori's (1997) 'median' flux (his eq. [3]). For Helium we use Engelmann et al. (1985: Fig. 11) 
(HEAO-A3), and the interstellar spectrum given by Seo et al. (1991: Fig. 11a). 

For the injection spectra of protons, we find 7 = 2.15 reproduces the observed spectra in the 
case of no reacceleration, and 7 = 2.25 with reacceleration. We find it is necessary to use slightly 
steeper (0.2 in the index) injection spectra for Helium nuclei (see Table [l]) in order to fit the 
observed spectra in the 1-100 GeV range of interest for positron production. The spectra fit up to 
about 100 GeV beyond which the Helium spectrum without reacceleration becomes too steep and 
the proton spectrum with reacceleration too flat; these deviations are of no consequence for the 
positron calculation. Although the nucleons are not the main subject of this study, we note that 
the reacceleration model reproduces slightly better the observed spectrum below 10 GeV/nucleon 
where the interstellar power law in momentum continues down to 2 GeV/nucleon before bending 
over ( |Seo et al. 1991 ). The adopted nucleon spectra successfully reproduce the 100-1000 MeV 



gamma-ray intensity from 7r°-decay ( [Strong et al. 19"96| , [Strong fc Moskalenko 1997a , [Strong 



Moskalenko, & Schonfelder 1997), so that we are not dependent only on the locally measured 
nucleon spectra for the positron calculation, but have evidence that it extends thoughout the 
Galaxy. Possible gamma-ray evidence for a flatter nucleon spectrum ( [Hunter et al. 19*9^ ) and the 
consequences for positrons are addressed below (Section 6). 

For the primary electrons, an injection index of 2.1 below 10 GeV, steepening to 2.4 above 10 
GeV, was found to reproduce the observed spectrum up to 30 GeV and is consistent with gamma- 
ray and synchrotron radiation studies ( Strong fc Moskalenko 1997a| ). At higher energies a further 



steepening is required but this is not of consequence here. The flux normalization was chosen to 
fit direct measurements and the positron ratio as described in section 5: I e (E) = 3.2 x 10~ 8 cm -2 
s- 1 sr- 1 MeV- 1 at 9 GeV. 



4. Production of secondary positrons and electrons 



The production of positrons in collisions of cosmic ray protons with protons of the interstellar 
medium has been discussed in detail in numerous studies (e.g., Stecker 1970, Orth &: Buffington] 
1976, Protheroe 1982, Dermer 1986a b, Murphy, Dermer, &; Ramaty 198"^ and reference therein). 
The muons created through decays of secondary pions and kaons are fully polarized, which results 
in electron/positron decay asymmetry, which in turn causes a difference in their production 



spectra. In the early papers ( Orth &: Buffington 1976 ) the muon decay asymmetry was treated 
identically for positrons and electrons (our paper: £ = — 1 in eq. flCl| ), which leads to almost a 
factor 2 lower positron production at high energies. Subsequent papers appear to have repeated 
this mistake (e.g., Protheroe 1982j ). The difference in /j + and \i~ decay asymmetry as applied 
to positron and electron production was noted by Dermer (1986a), but no attempts to improve 
the predicted cosmic-ray positron spectrum seem to have been made. Additionally, the energy 
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spectra of positrons produced in collisions of isotropic monoenergetic protons with protons at 
rest calculated by Murphy, Dermer, & Ramaty (1987), although in general agreement with our 
calculations, still differ in detail. In particular, they agree better with our calculations using an 
isotropic distribution of positrons in the muon rest system. Since positron/electron production 
and propagation in the Galaxy is among the hot topics of cosmic-ray physics we give explicitly our 
formulae for the positron/electron production spectrum (see Appendix). 



The energy spectra of positrons (£ = +1 in eq. []Cl]]), and electrons (£ = —1) from the decay 
of ir^ , mesons produced in collisions of isotropic monoenergetic protons with protons at rest 
are shown in Fig. ||| for several values of the kinetic energy of protons. For comparison we show 
the spectra of positrons calculated assuming an isotropic distribution in the muon rest system, 
£ = 0, and also using £ = — 1; the latter corresponds to the formula used by Orth & Buffington 
(1976). Remarkably, the latter one produces more positrons at maximum (around 30 MeV) as 
compared to the correct positron spectrum, while it gives half the positron yield at high energies. 
The curve corresponding to the isotropic distribution lies exactly in between the two others. 

The effect is also clearly seen when integrating over the spectrum of cosmic-ray protons (Fig. 
where we adopted for illustration the spectrum by Mori (1997), in units of cm~ 2 s _1 GeV : 
J p = 1.67p~ 2 - 7 (l + [2.5 GeV/c] 2 /Pp) _1/2 below 100 GeV, and J p = 6.65 x 10- 6 (e p /100 GeV)~ 2 ' 75 
above 100 GeV. The difference in the positron production spectra for the £ = +1 and £ = — 1 
cases reaches a factor 1.6 above ~ 0.2 GeV. It supports the conclusion made by Orth &: Buffington 
(1976) that the neglect of the muon decay asymmetry can result in a 25% error in the positron 
production spectrum, but curiously the inclusion of the correct kinematics leads to increase of the 
positron yield at high energies, just the reverse of their inference. The effect of kaon decay is of 
minor importance. We included only the most important channel, — > ^ + Vn, the others 
would contribute at the few per cent level. A comparison with the positron production rate from 
Protheroe (1982) illustrates the statement above: it generally agrees with our calculations for the 
case £ = — 1, also taking into account that his calculations include a contribution of He nuclei in 
cosmic rays and interstellar matter. Our curves are shown for pure pp- interactions, while inclusion 



of He would give a factor of ~ 1.4 increase (e.g., Dermer 1986b ). Some discrepancy could be 
connected with uncertainties in extrapolation of the inclusive cross section used by Protheroe at 
high energies and interaction dynamics. The uncertainty at low energies (hatched area) is due to 
uncertainty in the demodulation of the proton spectrum. 

Our calculations of the interstellar e^ spectra include the contribution to ir^ production from 
the channels p + He, a + H, and a + He. For collisions involving nuclei with atomic numbers A > 1, 
the corresponding cross section is multiplied by a factor (A\^ + A^ 8 — l) 2 (Orth & Buffington 



1976, Dermer 1986a), while the energy per nucleon is the kinematic variable. 
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5. Interstellar positron and electron spectra 

Fig. H shows the computed secondary positron and electron spectra for the cases without and 
with reacceleration. The primary electron spectrum and the total electron + positron spectrum 
are also shown. Above a few GeV where solar modulation is small the agreement with the absolute 
positron measurements is good within the large experimental errors. Unfortunately there are few 
absolute measurements of the positron spectrum, and we have to rely mainly on the positron 
fraction for comparison. First, in Fig. ^, we show the results plotted as the positron fraction using 
the total electron spectrum from Protheroe (1982). Next, Fig. [?] shows the positron fraction but 
using our computed interstellar primary electron spectrum. None of the computed positron spectra 
or positron fractions are seriously in conflict with the observations. Between 1 and 10 GeV the 
agreement is very good and the measured data rather precise. Above 10 GeV there appears to be 
an excess above the predicted ratio although the observational errors are larger. Below a few GeV 
solar modulation will shift the points to lower energy, but the implications for modulation models 
are beyond the scope of this paper. 

Reacceleration produces a more peaked positron spectrum with the peak positron fraction 
around 600 MeV compared to 400 MeV without reacceleration. Although the better fit to the 
positron fraction is given by the model without reacceleration (Fig. this depends critically on 
the measured electron spectrum which is not known to sufficient accuracy to allow a distinction 
between the models on this basis. While our non-reacceleration electron spectrum is consistent 



with Galactic gamma-rays from 1-1000 MeV and radio synchrotron data down to 38 MHz (Strong 



fe Moskalenko 1997a, Strong, Moskalenko, & Schonfelder 1997), the reacceleration spectrum may 



not be consistent with such data. It will in any case be possible to put critical constraints on 
electron re-acceleration using gamma-ray and radio data, and this will be investigated in future 
work ( [Strong fc Moskalenko 1997b ). Note that Protheroe (1982) uses a radio-based electron 



spectrum below 1 GeV, but this is increasingly uncertain at lower energies due to free-free 
absorption of radio emission. 

Our positron spectrum without reacceleration is steeper than that computed by Protheroe 
(1982); ours is a factor 2 higher at 1 GeV and is about equal to his at 10 GeV. A factor 2 and 
part of the steeper slope can be traced to a difference in the production function (see Fig. ||), and 
the remaining difference must be attributed to many details in the different propagation models 
(diffusion coefficient, halo size, gas density, radiation fields, 3-D versus 1-D models). Taking 
these reasons into account the calculations are consistent and it is gratifying that the difference 
is relatively small. However the new calculations are able to benefit from better knowledge of the 
nucleon spectrum and of cosmic-ray propagation. 

Above 10 GeV the predicted positron flux appears too low and the agreement is not better 
than for Protheroe (1982). This would be a hint for primary positrons or a harder interstellar 
nucleon spectrum (see Sect. 6) than observed locally and emphasizes the importance of accurate 
positron measurements in this range. At energies below 1 GeV where solar modulation is large 
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our interstellar positron spectrum should provide a good basis for modulation studies. 

As promised in the Introduction, we can use these results to turn the argument around and 
draw conclusions about the interstellar electron spectrum using the calculated secondary positrons 
and the measured positron fraction. In Fig. [7| we divide the computed positrons by the modelled 
primary electron spectrum plus the secondary e^\ the agreement with the measured positron 
fraction shows that either of the modelled primary electron spectra (Fig. |5|) are acceptable in the 
range 0.1-10 GeV. They are consistent with the direct measurements of the primary spectrum. 
The model without reacceleration is however lower by a factor of about 2 than that required 
to fit the COMPTEL and EGRET gamma-ray spectrum (see Strong fc Moskalenko 1997a|) ; an 



intermediate spectrum would nevertheless be consistent with both positrons and gamma rays 
within the uncertainties of each. The model with reacceleration may not be compatible with 
gamma-ray and radio data, as mentioned above. 



6. A harder interstellar nucleon spectrum ? 

The fit to the positron fraction above 10 GeV could be improved by an ad hoc steepening of 
the electron spectrum, but this would then disagree with the direct electron measurements which 
are rather reliable in this range. Another possibility which we find interesting is to adopt harder 
interstellar proton and Helium spectra. Fig. |8] shows the proton spectrum and positron fraction 
for an injection spectral index 7 = 2.0 (model 08-009, no reacceleration). The primary electron 
normalization for this case was adjusted to I e {E) = 3.44 x 10~ 8 cm~ 2 s~ x sr _1 MeV -1 at 9 GeV to 
obtain the best overall fit. This model reproduces the positron fraction well over the whole range; 
in particular the data above 10 GeV from the recent HEAT experiment ( [Barwick et al. 1997 ) are 



fitted (note that the higher positron fractions from earlier experiments can probably be attributed 
to the difficulty of distinguishing positrons from protons). The ambient proton spectral index after 
propagation is about 2.6 in this model, compared to the directly measured value of 2.75. This is 
especially interesting in view of the independent result from EGRET diffuse Galactic gamma-ray 
data ( Hunter et al. 1997| ) that the gamma-ray spectrum above 1 GeV is harder than expected 



for the locally measured nucleon spectrum (Mori "19971 , IGralewicz et al. 1997| ). Mori (1997) finds 



that an ambient spectral index of 2.41-2.55 is required to fit the gamma-ray spectrum, which is 
consistent with that required for the positrons. Thus our result can be taken as adding some 
support to the 'hard nucleon spectrum' interpretation of the gamma-ray results, and for the idea 
that the Galactic interstellar proton and Helium spectra are harder than those measured in the 
heliosphere. However since this situation is a priori unlikely in conventional models of cosmic-ray 
propagation and from cosmic-ray anisotropy arguments, it remains an interesting possibility 
deserving further investigation and independent tests. If correct it would explain the positron data 
without the requirement for a source of primary positrons. 
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7. Conclusions 

We have carried out a new computation of the secondary positron and primary electron 
spectra in a self-consistent model of propagation including nucleons. The model is more realistic 
than previous leaky-box type calculations and incorporates a wider range of astrophysical input. 
We have shown that the positron fraction is consistent with measurements up to 10 GeV, beyond 
which some excess is apparent. We have also shown that it is possible to reverse the normal 
argumentation to constrain the interstellar electron spectrum on the basis of the measured positron 
fraction and the computed positron spectrum. The resulting interstellar electron spectrum for 1-10 
GeV is consistent with that from direct measurements. A harder interstellar nucleon spectrum 
allows the positron fraction to be fitted also above 10 GeV and would also explain the high-energy 
diffuse Galactic gamma-ray spectrum; thus two independent lines of evidence point to a natter 
nucleon spectrum than that directly measured, so that this possibility has to be taken seriously 
and further consequences examined. At energies below 1 GeV, our calculated positron spectrum 
will be of interest for studies of cosmic-ray modulation in the heliosphere. 

More details, including the software and datasets, can be found on 



http:/ /www, gamma. mpe-garching.mpg.de/^aws/aws.htmi and ~zmos/ 'imos.html 



APPENDIX 

A. Spectra of secondaries from pp-collisions 

The production spectrum of secondary 7's, electrons and positrons can be obtained if one 
knows the distributions of pions F n (e 7T ,e p ) and kaons Fk(ek,£p) from a collision of a proton of 
energy e p , and the distribution of secondaries F s (e s ,e n) K) from the decay of a pion/kaon of energy 

j s {e s )=n H V / de p J p {s p ) (r]<ji(e p )} dei F s (e s ,ei)Fi(e h e p ) , (Al) 



,K 



where n# is the atomic hydrogen number density, J p (s p ) is the proton flux, and (jyov.A'^p)) is the 
inclusive cross section of pion /kaon production (a convenient parametrization for different channels 
is given by Badhwar, Stephens, & Golden 1977, and Dermer 1986a| ). The minimum proton energy 



e™ in that contributes to the production of a meson with energy £™ m , and the maximal energy 
of the produced meson e™ ax (e p ) can be easily derived from kinematical considerations, e.g. by 



equating y/s = [2m p (e p + m p )] 1 '' 2 = \Jm x + e* 2 — mf + e*, where s is the square of the total 
energy in the center-of-mass system (CMS), mx depends on the reaction channel, and e* is the 
CMS energy of the produced meson (connected with the laboratory system energy £j via the 
Lorentz transformation). The minimum meson energy that contributes to the production of a 
secondary with energy e s is given by e™ 11 ^) = e 7 + m 7r o/(4e 7 ) for 7-rays; £™™ K ± = fn- n ±^K ± f° r 
e e < E = ±7^7^(1 + /3p, and e™i n ^ ± (e e ) = ±m n , K (e e /E + E/e e ) if e e > E, where 7^, are the 
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muon Lorentz factor and speed in the meson rest system (see below). 



B. Pion production in pp-collisions 

We consider pion production in pp-collisions following a method developed by Dermer 
(1986ab), which combines isobaric ( [Stecker 1970| ) and scaling (Badhwar, Stephens, & Golden 
|1977 , [Stephens Badhwar 1981 ) models of the reaction. The isobaric model was shown to work 
well at low energies, while at high energies the relevant model is based on scaling arguments. In 
the transition region we join the models with a linear connection in the regime between 3 and 7 
GeV. 



B.l. Stecker 's model 

Assuming that the outgoing A-isobar of mass m A travels along the initial direction of the 
colliding protons in the CMS and decays isotropically, the distribution of the 7r's in the laboratory 
system (LS) is 

e p ;m A ) = 1 j * iJ[ 77r ; 7 + 7 ;(l - /3+#), 7 +Y w (l + 

' H[ln, 7at4(1 - /Mr), 7X7^(1 + ) , (Bl) 



7 a Pa 

where H[x; a, b] = 1 if a < x < b and = otherwise, e n is the pion energy, and e p is the LS energy 
of the colliding proton. The Lorentz factors of the forward (+) and backward (-) moving isobars 
in the LS are 7 ^ = 7c7 ^(l ± /3 c /3^), where 7c = y/s/2m p is the Lorentz factor of the CMS in the 
LS, 7^ = (s + rn\ — m p )/2^m / \ is the Lorentz factor of the isobar in the CMS. The pion Lorentz 
factor in the rest frame of the A-isobar is j' n = {m\ + — m p )/2mATnir- 

Integration over the isobar mass spectrum (Breit-Wigner distribution) yields the distribution 
of pions 

jp , \ r L -l (Vs-m p -m° A \ _! (m p + m,, - m A \ \ 1 

F ir (£ 1) £p) = r tan ~~ -tan 1 - 1 1 



r )\ 

where m A is the average mass of the A-isobar, and V is the width of the Breit-Wigner distribution 
(note that T here is a factor of 2 higher than the value usually given in particle data tables). 
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B.2. The scaling model 

The Lorentz invariant cross sections for charged and neutral pion production in pp-collisions 
inferred from experimental data at e p ^ 13.5 GeV are given by Badhwar, Stephens, & Golden 
(1977), and Stephens & Badhwar (1981) 

d 3 a 

£,,-=— = AG n (e P ){l ~ x n ) Q exp[-B P± /(l + 4m?/a)] , (B3) 



where 



G„±(e p ) = (l + 4m 2 p /s)- R , (B4) 
G n o(e p ) = (1 + 23e; 2 - 6 )(l - Am 2 p /s) R , (B5) 

Q = (Ci- C 2P± + C 3 p 2 ± )/Jl + 4ml/s , (B6) 



x n = ^xl + (4/s)(p 2 ± +ml), (B7) 
# 2m n ^/s -y c ~/ n (p n cos 6> - /3 C ) 



[(s — m 2 — m 2 ^) 2 — 4m 2 m 2 f ] 1 / 2 ' 



(B8) 



6* is the pion LS polar angle, A, B, Ci^,3, R are the positive constants, and rax depends on 
the reaction channel: rax = 2m p for reaction pp — > ir°X , mx = m p + m n for pp — > ir + X, 
mx = raj for pp — > 7r + (i, mx = 2m p + m^ for pp — > 7r _ X, mx = m p + m n for — > ET + X, and 
rax = 2ra p + m^ for pp — > ET~X. 

The LS energy distribution of pions can be obtained by integration over the LS polar angle 

F^e^^Ep) = r ^{- — / dcosO £^-3 — , (B9) 

where, provided —1 < cos# m i n < 1, 

cos e min = - J— f 7c£ . - s ~ m x + < \ (B10) 

and (^o"7r(^p)) sm is the inclusive cross section of pion production in the scaling model. 

The charged kaon production is of minor importance in comparison with pion production. 
The parameters of the Lorentz invariant cross section 

d 3 a 

£K-x— = ^(1 " *k) C exp(-Bp ± ) (Bll) 

are given by Badhwar, Stephens, & Golden (1977). The two main decay modes are —> fj,Un 
(63.5%) and — > ir ^ (21.2%), other modes, with three particles in the final state, contribute 
at a few per cent level ( particle Data Group 1990|) . 
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C. Pion decay 



The muons from the decay of charged pions are created fully polarized, which results in 
decay asymmetry. The secondary electron/positron distribution in the muon rest system assuming 
a massless electron/positron is given by ([Particle Data Group 199C| ) 



4(4 cos 9') 



8s't 



rrr 



3-4- 



m, 



£cos0' 1-4- 



m. 



(CI) 



where £ = ±1 for e^, £ = for the isotropic distribution, e' e , 6' are the electron/positron energy in 
the muon rest system and polar angle respectively. 



In the LS, the electron/positron distribution is 

ft(e e ,ej = f 1 dcos9fs(e' e ,cose')J f £ Al^\ , 
Jcos8 L s \e e ,cosOJ 



(C2) 



where cosOl = max{— 1, (1 — rriu/'ZzeYu)/ Pu}i j u , /?« are the muon Lorentz factor and speed in 
LS, J is the Jacobian 

e e . 1 



J 



e e ,cos6> / 4 7^(1 - (iy, cos (J)' 



(C3) 



and 



COS0' 



e e 7 M (l - P^cosd), 
cos 6 — P u 

1 — Pa COS 6* 



(C4) 



The LS distribution of fi from the 7r-decay is 



1 



(C5) 



2m ll <y 7T p 7T y fl P' u " " 

where £ M is the muon LS energy, 7^ = 7^7^(1 ± P^P'), 7' , ~ 0.2714 are the muon Lorentz 
factor and speed in the pion rest system, and 7^, p n are the pion Lorentz factor and speed in the 
LS. Therefore, the LS distribution of from the decay of pion with Lorentz factor 7^ is given by 

1 

F$ (£ e , e w ) — — 

m^P^P'^ 

~ < £ e < m M /2 7 +(l + P+) 

X 5 ( 71 ) - X c ( 7 -) + y c ( 7 +) - y f ( 7 i),m M /2 7 +(l + /?+) < £ e < m„/27-(l + /?") 



X \ 1 

[y c ( 7 +)-n(7i), 

where 71 = ^ + jg, and 



m^/2 7 -(l +/?-) < £ e < m u /2 lu (1 - /£) 
m M /2 7 -(l -/?")< £ e < m M /2 7 +(l - /?+) 



(C6) 



dcos f'^(s' e , cos 0') J 



e' , cos 0' 



£ e , COS 



J_fl 

/3 M \ 2e e 7/j 



(icOS fUt'e, COS 0') J 



£g, COS 0' 
£ e , COS 9 



(C7) 
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After the integration, one can obtain 



*±(7) = h ( — J | — [-32 7 3 (1 ± P) + 7(24 ± 32/3)] + 7 2 (27 ± 9/3) =f 91n[ 7 (l + /3)] 1 , (C8) 



9 \ m 



m 



*o(7) 



m. 



(C9) 



12 



2 



1.0 In ( } _ 647(1 -/9) 



7-1 



+ 



m. 



48 7 2 (1 -/3) + 24 In 



/3 



1 + /3 



21n( 7 /3) + 101n[ 7 (l + /3)] 



(CIO) 



y_( 7 ): 



36 I U 



-512 7 3 (1 -0)+ 7 (576 - 320/3) - 48 In 



7+1 



2 r 



+ 



288 7 2 (1 - /3) - 72 In 



P 



1 + /3 



7-1, 

+ 61n( 7 /3) + 301n[ 7 (l + /3)] 



(Cll) 



n(7) 



1 - 
18 ] I m, 



-128 7 3 (1 -13) + 7 (96-32/3)] 



+ 



m. 



108 7 2 (l-/3) + 151n[ 7 (l + /?)U, (C12) 



At large Lorentz factors, 7 ^ 20, it is necessary to use the series expansions. 
The distribution of 7 -rays from ir° decay is given by 

2 1 1 

F 7 (e 7 , e^o) = — , -m 7r77r (l - /3 7r ) < e 1 < -m^(l + f3 w ), 

Pit & 



(C13) 



where the factor two accounts for two photons per decay. 



D. Kaon decay 

The two main kaon decay modes are — »■ iiv^ (63.5%) and K ± — »■ vr ^ (21.2%), and 
thus they also contribute to the secondary e ± spectrum. The first mode is similar to 7r =l= decay, 
and all formulas presented in previous section are valid if one replaces index ir with K and takes 
appropriate values for 7 ' 0' . 
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The charged pion distribution from the two-pion decay of the kaon is 

j r w (e WJ sk) = r \ — j-sp, <e n < m w ^, (Dl) 

where is the pion LS energy, 7^ = 7^7^.(1 ± "/ ±, ft ± « 0.8267 are the charged pion 

Lorentz factor and speed in the kaon rest system, and "fx, Pk are the kaon Lorentz factor and 
speed in the LS. Therefore the LS distribution of pions from the reaction pp — > K^X is given by 

F n (e n ,e p )= 7rm . jK 3 . {ri(TK{e p ))~ l [ de K [ dcosO ( e K 4r^- ) , (D2) 

where cos # m i n is given by eq. ( BlOf ) with the replacement of the index tt with K, and e™ n (e 7r ) can 
be obtained from the equation e n = 7^-(e™ n + /%.p™ n ). 
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Table 1. Parameters of models. 

















protons 




Helium 




Model 


z h , kpc 


-Do, cm 2 /s 


Po, MV/c 


S 


va, km/s 


7 


Po a 


Io h 


7 


Po a 


/o b 


08-005 


3 


2.0 x 10 28 


3.0 x 10 3 


0.60 





2.15 


10 4 


3 x 10~ 6 


2.35 


4 x 10 4 


4 x 10~ 8 


08-006 


3 


4.2 x 10 28 


3.0 x 10 3 


0.33 


20 


2.25 


10 4 


3 x 10~ 6 


2.45 


4 x 10 4 


4 x 10~ 8 


08-009 


3 


2.0 x 10 28 


3.0 x 10 3 


0.60 





2.00 


10 4 


3 x 10~ 6 


2.00 


4 x 10 4 


4 x 10 -8 



a in units MeV c 1 nucleus 1 

b in units cm~ 2 s _1 sr _1 (MeV/c/nucleus) _1 

Note. — Model 08-005 and 08-009 are without reacceleration, model 08-006 is with reacceleration. Spatial 
diffusion coefficient for no reacceleration: /3Do below rigidity po, (3D (p/po) s above rigidity po- Spatial diffusion 
coefficient with reacceleration: /3Do(p/po) S for all rigidities. Nucleon injection density spectrum: dq(p)/dp oc p 1 . 
Nucleon spectra normalized to flux Jo at momentum po- 
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Fig. 1. — B/C ratio. Thick line: model with no reacceleration, nucleon injection spectrum index 
2.35. Thin line: model with reacceleration, nucleon injection spectrum index 2.25. Dashed lines: 
modulated to 500 MV. Data: Voyager, Webber et al. (1996). 
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Fig. 2. — Left panel: spectra of protons. Thick line: model with no reacceleration, injection 
spectrum index 2.15. Thin line: model with reacceleration, injection spectrum index 2.25. Dashed 
line: Seo et al. (1991) 'interstellar'; Dash-dotted line: Mori (1997) median spectrum. Right panel: 
spectra of Helium. Thick line: model with no reacceleration, injection spectrum index 2.35. Thin 
line: model with reacceleration, injection spectrum index 2.45. Dashed line: Seo et al. (1991) 
'interstellar'; Dash-dotted line: Engelmann et al. (1985). 
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Fig. 3. — Left panel: solid lines show the energy spectra of positrons from the decay of 7r + , K + 
produced in collisions of isotropic monoenergetic protons with protons at rest for various proton 
kinetic energies (from bottom to top): e v - m p = 0.316, 0.383, 0.464, 0.562, 0.681, 1.0, 1.78, 3.16, 
10.0, 100.0 GeV (£ = +1, eq. fCljl ). Dashed lines show the spectra calculated assuming an isotropic 
distribution of e + in the muon rest system (£ = 0) for some proton kinetic energies only. Dotted 
curves show the spectra calculated for £ = — 1, corresponding to the distribution used by Orth & 
Buffington 1976, their eq. [D9]. Right panel: energy spectra of electrons from the decay of ir~ , K~ 
produced in collisions of isotropic monoenergetic protons with protons at rest for several proton 
kinetic energies (from bottom to top): 1.0, 1.21, 1.47, 2.15, 3.16, 4.64, 10.0, 21.5, 46.4, 100.0 GeV 

(£ = -!)• 
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Fig. 4. — Illustration of the effect of the different positron distributions in the muon rest system. 
The thick solid lines show the production spectrum of secondary positrons and electrons per 
hydrogen atom for the cosmic-ray proton spectrum given by Mori (1997). The positron spectrum 
with no kaon contribution is shown by the dash-dotted line. The dashed line shows the spectrum 
calculated assuming an isotropic distribution of e+ in the muon rest system (£ = 0), the dotted 
line shows the spectrum calculated for £ = — 1. The thin solid line with hatched regions shows 
the positron production rate by cosmic rays in interstellar medium (including contribution of He 
nuclei) from Protheroe (1982). 
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Fig. 5. — Spectra of secondary positrons and electrons, and of primary electrons. Left panel: model 
with no reacceleration (08-005). Electron injection spectrum index 2.1, 2.4 above and below 10 
GeV respectively. Model: upper curves: primary electrons and primary+secondary electrons and 
positrons. Lower curves: secondary positrons, electrons. Lower dashed-dot line: Protheroe (1982) 
leaky-box prediction. Data: electrons: Buffington, Orth, & Smoot (1975), Golden et al. (1984), 
Golden et al. (1994), Taira et al. (1993), Ulysses ( Ferrando et al. 1996[ ), upper dashed-dot line: 
Protheroe (1982); positrons: Fanselow et al. (1969), Buffington, Orth, & Smoot (1975), Golden et 
al. (1987), Golden et al. (1994). Right panel: same, model with reacceleration (08-006). 
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Fig. 6. — Left panel: Positron fraction for model with no reacceleration. The positron spectrum is 
divided by the electron spectrum used by Protheroe (1982). Dot-dashed line: positron fraction from 
Protheroe (1982), leaky-box (dashed-dot), diffusion (dashed). The collection of the experimental 
data is taken from Barwick et al. (1997). Right panel: Same positron fraction for model with 
reacceleration. 




Fig. 7. — Left panel: Positron fraction for model with no reacceleration. The positron spectrum 
is divided by the electron spectrum computed in the propagation model. Data and other curves as 
Fig. H Right panel: Same positron fraction for model with reacceleration. 
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Fig. 8. — Left panel: spectrum of protons for a 'flat' injection spectrum. Thick line: model with 
no reacceleration (08-009), injection spectrum index 2.0. Other spectra as in Fig. Right panel: 
positron fraction for this nucleon spectrum. Data as Fig. ||. 



